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Abstract 

We show that sound waves scattered from a hydrodynamic vortex may be amplified. Such 
superradiant scattering follows from the physical analogy between spinning black holes and hydro- 
dynamic vortices. However a sonic horizon analogous to the black hole event horizon does not exist 
unless the vortex possesses a central drain, which is challenging to produce experimentally. In the 
astrophysical domain, superradiance can occur even in the absence of an event horizon: we show 
that in the hydrodynamic analogue, a drain is not required and a vortex scatters sound superradi- 
antly. Possible experimental realization in dilute gas Bose-Einstein condensates is discussed. 

PACS numbers: 03.75.Nt, 03.75.Kk, 04.62. +v 
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Hawking's prediction of radiation from black holes is one of the most surprising outcomes 

fl 

of the quest to unify quantum mechanics and general relativity It is problematic that 
such an important theoretical prediction has no immediate prospect of direct experimental 
investigation. Consequently, Unruh's observation of an analogy between wave scattering 
from black holes, and sound wave scattering from sonic horizons in fluids, has intrigued 
both the astrophysical and condensed matter communities □ ODD 

The generation and detection of the sonic analogue of Hawking radiation faces many 
difficulties, and a detailed experimental scheme has yet to be presented. Consequently, we 
propose the initial investigation of the simpler, but closely related, superradiant amplifi- 
cation of scattered waves by an ergoregion f|. Quantum mechanically, Hawking radiation 
and superradiance may be regarded as spontaneous and stimulated radiation counterparts. 
However, stimulated radiation is easier to observe because of its much greater strength, and 
because of the control provided by the stimulating pulse. Furthermore, an analogue ergore- 
gion is realized by a vortex, which is experimentally simpler to achieve than a horizon. This 
superradiance is related to the Zeldovich-Starobinsky effect, which has been discussed in the 
context of superfluid helium by Volovik Q|. We provide quantitative details and numerical 
simulations to argue that this superradiance might be observed in dilute gas Bose-Einstein 
condensates (BECs) Q,[g, 10, 11 1. 



According to general relativity, rotating black holes have an ergoregion: a region in 
which no physical object can remain at rest with respect to an inertial observer at infinity. 
Low frequency waves are scattered from the ergoregion with increased amplitude; that is 
superradiantly jl^. In fluids, regions where the fluid speed exceeds the speed of sound are 
analogue ergoregions 0| . A draining vortex with a central sink may also possess an analogue 
event horizon, and hence is an analogue of a rotating black hole. Basak and Majumdar [l 3l | 
showed that sonic superradiance occurs for a draining vortex, but assumed that the density 
of the fluid, and the speed of sound, were constant. 

For a BEC a sink might be realized by outcoupling the condensate to form an atom laser 
beam Q]. However there are difficulties with strong outcoupling and it will be simpler 
to create a vortex without a sink, for example by stirring with laser beams 16] or magnetic 
fields Then the simpler velocity function makes it feasible to consider the variable 

density profiles which are unavoidably associated with vortices in BECs. 

The wave equation for sound waves, that is linear perturbations of the velocity potential 
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of a barotropic, inviscid, irrotational fluid, may be written 



(1) 



where the coefficients g^ u are functions of the unperturbed fluid density and velocity, and 
q is their inverse determinant. This has the same form as the Klein-Gordon equation for a 
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massless scalar field in curved spacetime with metric g^ v |18l| . Provided the quantized fluid 
has suitable commutation relations, the analogy extends to the quantum field jlflj ]. 

All this suggests the possibility of using fluid flows to mimic general relativistic spacetime 
metrics. The analogy requires the length scales of the problem to be large enough that the 
atomicity of the fluid can be neglected. In BECs, furthermore, the hvdrodynamic equation 
fll} is an approximation neglecting the "quantum pressure" term 20(. This is valid when 
any spatial variations in the BEC occur over length scales longer than the healing length of 
the condensate (21]. In our case, this must be smaller than both the sound wavelength and 
the spatial scale of the vortex. 

Let the unperturbed fluid's speed of sound be c, velocity be v, and density be p. In 
cylindrical polar coordinates a vortex centred on the origin, with v r = v z = 0, has wave 
equation coefficients 
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and g = —p 4 r 2 /c 2 . For irrotational flow, vg — a/r, for some constant a. The density and 
speed of sound depend only on the radial coordinate r, and approach asymptotic values p^ 
and Cqo for r — > oo. 

We consider cylindrical wave solutions to the wave equation of the form <p(t, r, 8, z) = 
if)(t, r)e~ ime , with angular wavenumber m. Assuming that the square of the speed of sound 
is proportional to the density, as is the case for a BEC j^, the density may be eliminated 
from the wave equation, which then becomes 



d 2 %b mvg dip 
— 2i- 



r dr V dr 
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For single frequency waves of the form tp(t, r) 
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(C 2 -V 2 g)lP = 0. 



(3) 
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G{r*)e lu)t , and defining a "tortoise 



coordinate" r* by dr/dr* = c 2 , with r* — > r for r — > oo, where c = c/cqo, we find 



+ |v - VWGM = o , (4) 
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where we have introduced an effective potential V e ff to emphasise the similarity to the time- 
independent Schrodinger equation. 

In order to analyse superradiance we consider the two limiting cases where the tortoise 
coordinate r* approaches ±oo. For r large, r* « r, and eliminating all terms except those 
of highest order in r yields the asymptotic form 



d 2 G(r*) 
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+ — G(r*) « , (6) 



P( r ) = P™^~T, 7T12 ■ ( 8 ) 



dr* 2 c 2 ^ 
with general solution 

G(r*) = Ae i{uj/Cac)r * + B e - i{uj/c °° )r * , (7) 

where A and -B are constant amplitudes of incoming and outgoing waves respectively. 

The analysis of the other limit, where r* — > — oo, requires that the density profile of the 
vortex be specified. In order to present analytical calculations we use the density profile 

[jr -r )/a] 2 
• 2 + [( r -r )/af 

This is similar to the charge I = 1 vortex density profile for a BEC [20]. but with the 
scale length given by the free parameter a, rather than by the healing length x- m the 
Thomas-Fermi limit of dominant atom interactions, the length scale for a charge I vortex, 
with angular momentum lh per particle, is l\ [22] . Although a vortex with I > 1 is unstable 
to decay into I single charge vortices, it may be stabilised by a pinning optical potential J3] . 
Since the density profile Eq.(JHJ) was chosen primarily for analytical convenience, we do not 
expect it to be exactly achievable experimentally, although something close to it should be. 
For example, the dipole potential generated by far detuned light has been used to engineer 



density profiles in BECs |24|. The other free parameter is the radius of the zero density 
core ?~o, which might also be engineered by applied potentials. Since we have assumed that 
the square of the speed of sound is proportional to the density, c 2 = p/poo, and the tortoise 
coordinate is 

1 , 2a 2 

r" = I dr = r . (9) 

cr r — r 



For r m ro, r* « — 2cr 2 /(r — ro) and c 2 ~ c^r — ro) 2 /(2<7 2 ). Using these approximations 
and retaining only terms of lowest order in r — ro, we find that asymptotically 

^G(r-) , W (10) 



dr* 2 c 2 



Vt = uj — ma/r . (11) 

This is an Euler-type equation and has solutions of the form G(r*) = \r*\ 5 , where 

5= ^±7, 7=^ign(fi)(l-8a 2 fi 2 /c 2 ) 1 / 2 . (12) 

The general solution to the asymptotic form of the wave equation is then 

G(r*) = \r*\ 1/2 (C\r*p + D\r*\^) , (13) 

where C and D are constants. According to Eq. (jl2j) . 7 is either real or purely imaginary. 
For 7 imaginary this asymptotic solution is oscillatory and may be written in the form 

G(r*) = |r*| 1/2 (Ce 711117 '* 1 + De~~< lnlr * 1 ) . (14) 

By making |r*| the argument of a logarithm we have anticipated that r* will be made 
dimensionless. This last form allows us to identify the two linearly independent solutions 
as outgoing and ingoing waves, with amplitudes C and D, respectively. Note that in this 
asymptotic limit r* < and hence In |r*| increases with decreasing r*. Note also, that it is 
the sign of the group velocity, not of the phase velocity, which determines the identification 
of outgoing and ingoing waves 25], and the sign of the group velocity is independent of the 
sign of Q. 

Let the Wronskians (of the solution and its complex conjugate) corresponding to the 
asymptotic solutions Eqs. (0[T3J) be denoted W(oo) and W(— 00). Then 

ty(oo) = — (\B\ 2 -\A\ 2 ) . (15) 

Coo 

If 7 is imaginary, 

W{-oo) = 2 7 (|C| 2 - ID] 2 ) . (16) 

n 

Following Vilenkin [25], consider a solution to the wave equation for G(r*) representing a 
wave originating at r* = +00 and having the asymptotic forms, for imaginary 7, 

e i(u/coo)r* 1 ft e -i(u/ccc)r* r * } qq 

G ^=< , ll/2 , ,*, ' ' • ^ 

T r * l/2 e - 7 ln|r * _^ _ oc 



Then by Abel's Theorem, the Wronskian of this solution is constant, since there is no first 
derivative term in the differential equation for G(r*), Eq. (0J). Thus, setting A = 1 and 
B = R in Eq. (TSj), and C = and D = T in Eq. (JED, and equating W{oo) = W(-oo), 

^(| J R| 2 -1) = -2 7 |T| 2 , (18) 

Coo 

and therefore, 

|i2| 2 = l-sign(fi)| 7 ||T| 2 Coo /cu. (19) 

In particular, if Q is negative, and T^O, then \R\ > 1, and superradiance occurs, since then, 
in the asymptotic region far from the vortex, the amplitude of the reflected wave exceeds 
that of the incident wave. This superradiance inequality Q < expands to w < majr\ , 
which requires ma > 0, so that the cylindrical waves propagate in the same direction as the 
vortex flow. 

In contrast to the draining vortex examined by Basak and Majumdar [y| arid the rotating 
black hole [^J, here there are two independent necessary conditions for superradiance: 7 
must be purely imaginary and Q must be negative. Combining these conditions gives a 
frequency-independent condition for superradiance. For since u > 0, then Q > — ma/r^, 
and superradiance requires, > Q > —ma/rQ, and therefore Q 2 < m 2 a 2 /rQ. If 7 is to 
be imaginary, we must have 1 < 8a 2 Q 2 /c 2 ^, and since ma > 0, the frequency-independent 
necessary condition for superradiance is ama > CooTq/^a/^). 

For a charge I vortex in a BEC we can substitute expressions |2Q| for the vortex velocity 
constant a = IS/m atom , where m atom is the mass of an atom, and for in terms of the 
healing length x, Cqo = ^/(-\/2x m atom), to obtain x > r^/ (Aalm). This is always possible 
to fulfill, at least in principle, since the healing length is given by x = (87ma)~ 1//2 , where n 
is the atomic number density and a is the s-wave scattering length, and both n and a are 
under experimental control. 

To determine the magnitude of the superradiant amplification, we numerically solved the 
wave equation Q using the XMDS package j^J- We used an initial Gaussian wavepacket 
of the form ip(0,r) = A(r), with frequency Uq. Specifically, 



— g-(r-'"tait) /w gi^WCoo 

^(0, r) = (ilu - 2 Coo (r - r init )/w 2 ) A{r) . (20) 
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FIG. 1: Wavepacket propagation from numerical solution of Eq. @. The real part of the wave 
packet is plotted versus the dimensionless modified radial tortoise coordinate r*, defined in the 
text, with e = 800/801 so that f* = —800 at r = Tq. The initial wavepacket, at t = 0, was given 
by Ea. pnj) with 

Hnit — 200 and w — 10. (a) t — 114: the wavepacket is propagating towards the 
vortex, (b) t = 330: the wavepacket has just split into reflected and transmitted parts, (c) t = 504: 
the reflected wavepacket. Dimensionless parameters: m = 1, a = 2, Cqo = 1, a = 20. The spatial 
grid had 8192 points over r* G [—800, 400], and the time steps were 0.015 dimensionless time units. 

We found that superradiance occured for a wide range of parameters, as long as the relevant 
inequalities were fulfilled. Fig. ^ shows a particularly strong and clean example of the 
scattering of an m = 1 wavepacket. Comparing the Fourier components of the real parts of 
the incident and reflected wavepackets in the asymptotic region, we found that the dominant 
Fourier power in the reflected wavepacket is approximately doubled. 

The parameters are given in the caption, and are representative of conditions for a trapped 
dilute gas BEC. They are made dimensionless, indicated by a tilde, by measuring distance 
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and time in units of r$ and u;^ 1 respectively. For example, choosing tq — 1 /iin and Uq 1 = 
1CT 3 s, the parameters of Fig. 1 correspond to Coo = 10~ 3 ms _1 and a = 20 /im. For the case 
of a Rb 87 BEC, m atom = 1.4 x 10~ 25 kg, and the asymptotic healing length x ~ 0.5, and 
a = 2 corresponds to a vortex charge of about / = 3. Note that the real part of the wave is 
plotted against the modified tortoise coordinate f* — f — 2a 2 / (f — 1 + e). This magnifies the 
scale as r — > r (f — * 1) , while limiting the lower bound of r* to 1 — 2a 2 /e. The ergoregion 
boundary, or static limit, is at f pa 8 or r* ~ —92. 

The local healing length determines the validity of the hydrodynamic approximation, 
and is larger than the asymptotic value by the factor a/ Pool p( r )- Because of this the BEC 
interpretation fails for Fig. 1(b). However, we have found a tradeoff between the strength 
of the superradiance and the validity of the BEC hydrodynamic approximation. With in- 
creasing m, a, and sound wavelength, reflection occurs at larger values of r, and hence for 
smaller local healing lengths, while the superradiance decreases in strength. We define the 
reflection point to be at the maximum of the effective potential Eq. ((S]). Compared to the 
case of Fig. 1, a transmitted wavelength comparable to the local healing length at the re- 
flection point, accompanied by a 70% reduction in superradiance power, occurs for an initial 
wavelength an order of magnitude longer, and a a factor of five larger. 

Numerical simulations of BECs under the proposed experimental conditions, and without 
the hydrodynamic approximation, could explore the effect of its breakdown as the wavelength 
of the transmitted wave decreases due to its inward propagation. This is related to the 
"trans-Planckian" problem of Hawking radiation [27 1 . 

In summary, our analytical and numerical work has shown superradiant sound scattering 
from hydrodynamic vortices. We have also argued for its experimental realization in dilute 
gas Bose-Einstein condensates. This might provide a useful step towards the ultimate goal 
of observing the sonic analogue of Hawking radiation. 

ACQAO is an Australian Research Council Centre of Excellence. This research was 
supported by the Australian Partnership for Advanced Computing. 
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